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We develop a free-carrier theory of the optical absorption of light carrying orbital angular mo- 
mentum (twisted light) by bulk semiconductors. We obtain the optical transition matrix elements 
for Bessel-mode twisted light and use them to calculate the wave function of photo-excited electrons 
to first-order in the vector potential of the laser. The associated net electric currents of first and 
second-order on the field are obtained. It is shown that the magnetic field produced at the center 
of the beam for the £ = 1 mode is of the order of a millitesla, and could therefore be detected 
experimentally using, for example, the technique of time-resolved Faraday rotation. 
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PACS numbers: 78.20.Bh,78.20.Ls,78.40.Fy,42.50.Tx 

It is well known from classical clcctromagnctism that 
light can carry spin and orbital angular momentum. 
While the former was detected for the first time in the 
1930s [l[, the latter became widely available for exper- 
imental study only recently after the work of Allen et 
al. [3] In a seminal paper, those authors showed that 
light carrying an integer amount of orbital angular mo- 
mentum {hl^ with I an integer) may be generated in the 
laboratory using conventional laser beams. Since then, 
research on the subject of light carrying orbital angular 
momentum (0AM), or twisted light (TL) 0, 0, Q has 
spanned a large number of areas, namely, generation of 
beams [J], interaction with mesoscopic particles (optical 
tweezers) [y, [Zl, |8|, entanglement with spins for poten- 
tial applications in quantum information processing [9|, 
interaction with atoms and molecules [lOl . [ll|, cavity- 



and interaction with Bose-Einstein condcn- 
lj|. Nevertheless, the interaction with solid 



QED 

sates 

state systems, although potentially important for tech- 
nological applications, has not been explored so far. In 
this Letter, we present the first theoretical predictions 
about the interaction of TL with bulk semiconductors. 
We consider band-to-band transitions, i.e. optical transi- 
tions with light frequencies above the bandgap, so that 
free carriers rather than excitons are produced. We show 
that there is a transfer of 0AM between the light and the 
photoexcited electrons so that a net electric current ini- 
tially confined to the beam area appears. The magnetic 
field induced by these photocurrcnts is estimated. 



A beam of TL presents an azimuthal phase depen- 
dence — helical wavefront — responsible for the 0AM, 
and a radial dependence of the Laguerre-Gaussian{hG) 
or Bessel mode type. We will focus on the Bessel modes, 
but our results are applicable to the LG-mode with slight 
changes. The vector potential in the Goulomb gauge with 



cylindrical coordinates {r|| , cj), z} is [la] 

Aq,±(r,<) = Aoe*(?^^-*) [e±Ji(<Z||r||)e*'*T 



~,n 



iz^Ji±i{qn)e<^^'^^ 



+C.C., (1) 



with polarization vectors e± = x ± iy, Bessel functions 
J;, and parameters qu '^ Qz- 

Semiconductors are solids which at zero temperature 
have the highest occupied (valence) and lowest empty 
(conduction) energy bands separated by a gap Eg. In 
this respect, they are closer to insulators than to met- 
als. However, in typical semiconductors i?g ~ 1 eV, 
making possible the transitions between the valence and 
conduction bands by optical excitation. In a crystalline 
semiconductor, electrons in the valence and conduction 
bands, denoted by A = {v,c}, occupy the Bloch states 
</'Ak(r) = L~'^/^e*''''uAk(r), with uxk{r) a cell-periodic 
function (lattice constant a) , /i k the crystal momentum, 
and L the linear size of the semiconductor. This sim- 
plified two-band model grasps the main features of a 
semiconductor [l6[, and can be considered a good model 
of a real bulk system under the following two conditions: 
i) an applied strain splits the heavy-hole and light-hole 
valence bands; ii) only circularly polarized light is con- 
sidered, producing optical transitions of only one electron 
spin type between valence and conduction bands. 

We now proceed to develop an analytical description 
of the interband transitions in a two-band model of a 
bulk semiconductor induced by twisted light. Our analy- 
sis tackles the coherent optical excitation to conduction- 
band states, i.e. we consider the case Hu > Eg. In 
this regime the creation of excitons is negligible and one 
only has to take into account the free carriers trans- 
ferred by optical excitation from the valence to the con- 
duction band [l6|. The lowest-order contribution to 
the light-matter interaction, using the minimal-coupling 
Hamiltonian, is Hj = —{Q/m) p • Aqic,{r) with Aq;cr(r) 
the transverse or xy part of Aq/±(r,f) [Eq. ^], and 
{p, 771, Q} the momentum operator, mass, and charge of 



the particle involved. The complete semiclassical Hamil- 
tonian — operators for electrons and classical variables 
for the light field — consists then of two terms, the bare 
electron energy and the interaction Hj 



Ak 

^ (A'k'lp- Aqi„(r)|Ak}a^,^,aAk: 



A'k'Ak 



with ali,,/axu the creation/annihilation operators for 
electrons in a Bloch state of band A = v,c 
and quasi- momentum ?ik, and Ex^ its bare energy. 
The action of p = —ih'V onto the wave func- 
tion 95Ak(r) = (r |Ak) yields two terms: piy9Ak(r) = 
fik(pAk(r) — ihL^'^^^e^^''^\7u\k{r). For interband tran- 
sitions - our interest - the first term is negligible; then 

ff I th Q f i(k'_k).r 

-n/|A'k'Ak — T^— / «re 



^A'k'(r)[Aq/<.(r)-V]7.Ak(r). 



This expression can be simplified by using a standard 
procedure in solid state physics. The integration over 
the complete system is replaced by an integration over a 
unit cell and a sum over all cells 



Hr\\iu' 



/lA'k'Ak 



raN 



^^-.(k'-k).R.A^^^(j^^ 

C 

— / dru^,k-(r)[-iftV]uAk(r), 



where N is the total number of unit cells in the crystal. 
The integrand is split thanks to the facts that, in the 
length of the unit cell a, Aqicri^) and exp [— i(k' — k) • r] 
are almost constant and u(r) is periodic. The integral 
is the momentum matrix element PA'k'Ak- The sum is 
handled using the Jacobi- Anger identity, which allows us 
to separate the angular from the radial dependences in 
g-t(k-k)Rc^ and then applying the orthogonality rela- 
tion for Bessel functions [l7|. The final expression for 
the interaction Hamiltonian is 



Hr 



'(-*) 
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e'^' (e<j ■ Pck'»k) a^k'^i-k + ^-C- 



(2) 



with K = k — k having cylindrical coordinates 
{k\\,0,K;,}. 

In order to understand the effect that a twisted-light 
beam has on the ground state of a semiconductor, we 
focus our attention on the positive part Hj of the in- 
teraction Hamiltonian - first term of Eq. ^. The action 
of Hj on the full ground state of the A^-electron semi- 
conductor can be shown to yield an eigenvector with ex- 
pectation value of the orbital angular momentum equal 




FIG. 1: Schematic representation of interband optical exci- 
tation with twisted light. A twisted-light beam excites a 
valence-band electron into a superposition of states of the 
conduction band in a "cone" -like fashion. The laser parame- 
ter gii equals the cone radius at the conduction band. 



to hi. Here we present a simplified version of the theory 
where we show the action of the Hamiltonian on a single- 
particle valence-band electron state. The analysis is re- 
stricted to optical excitations that couple electron wave- 
numbers near k = 0, a typical situation in semiconduc- 
tor optics. We define ^ = ^(^i)' {Q Ao/m) {eg- ■ PcOuo), 
transform variables {k, k'} -^ {k, k}, fix the value of 
k and take the continuum limit for k. In the coordi- 
nate representation and rotating frame, the action of the 
Hamiltonian on the valence band single-particle initial 
state gives 



ijj+V.k(r) ^ ^ [ dK 
Ztt J 



n 



5{kz - gz)e'*Vck+K(r) , 



which makes evident that a "cone" -like transition oc- 
curs producing a linear superposition of conduction-band 
states with varying phases, as depicted in Fig.[T] Exploit- 
ing the delta functions, using iy5Ak+K(r) — e*'^''(y5Ak(r), 
and applying the Jacobi- Anger identity 



i/r V.k(r) = i [i'e'«^V,(<Z||r||)e"^'] ^ek(r) 
= C./(r)¥'ck(r)- 



(3) 



oi4>l 



The dimcnsionless function /(r) = i e'''^' ^ Ji{q\\r\\)e'^ 
contains all the relevant information that distinguishes 
the action of a twisted beam from its counterpart, the 
plane-wave. Because of the presence of /(r), the RHS 
of Eq. ^ is not a conduction-band eigenstate. For the 
expectation value of the z-componcnt of the orbital an- 
gular momentum operator L^ = —ihd^ in the state of 
Eq. ([3]) we obtain < L, >= hl^ since < L^ >= hi for 
the state e**^' and < L^ >= for the states Uck(r) and 
gtkr^ Thus, for short times, the evolution leads to a well- 
defined transfer of orbital angular momentum hi from 
the light beam to the electron. 



This net transfer of orbital angular momentum to the 
photo-excited electrons is expected to result in electric 
currents and associated magnetic fields. The latter would 
be an experimentally detectable signature of the type of 
optical excitation described here, and may also lead to 
opto-electronic applications. An estimate of the total 
current/magnetic field produced by all electrons can be 
obtained by calculating the total number of electrons ex- 
cited by the field, then determining the current/magnetic 
field of a single electron, and, finally, multiplying the 
single-electron current/magnetic field by the number of 
photo-excited electrons. Wc note that this is a single- 
particle calculation which docs not take into account the 
electron-electron Coulomb interaction. This is justified in 
the present study by the fact that the main correction in- 
troduced by the Coulomb interaction would be the renor- 
malization of the single-particle energies [l8|, ll9| . This 
mean-field effect would be small since we are working in 
the regime of low excitation (the density of photogen- 
erated electrons and holes is small), and, furthermore, 
by its nature, it does not affect the qualitative picture 
drawn in this work. Beyond the mean-field approxima- 
tion, the Coulomb interaction introduces complex and 
possibly interesting scattering effects ISl, l20| , which arc 
beyond the scope of our present analysis, and are left for 
future study. Again, however, these effects are expected 
to be small in the low-excitation regime. 
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FIG. 2: First-order electric current for OAM ? = 1 to 4. The 
center of each plot coincides with the center of the twisted- 
light laser beam. 

The number of electrons excited in a volume dV and 
a time interval t is dNex = a I(r)t/huj, with a the ab- 
sorption coefficient, /(r) = Ioe~"^ is the intensity of 
the light field, which is related to the vector potential 



by 2/o = veuj"^ A^, where v is the speed of light in the 
medium. Then, the total number of electrons that are 
excited in V = L^ is 

with f(quL) a function containing J;_i(x), Ji{x), 
Ji+i{x), and e^"^. 

In order to calculate the electric current, we first de- 
termine the single photo-excited electron state. For short 
times such that t{u! + Eg/h) < 1 the evolution of the state 
under the action of the laser field is, in the Schrodinger 
picture, 

V'(r,i) = U{t)ip,i,{r) 

o, ^,k(r) - e-'^^*/'' {it/h)U{r)Vc^.{r) , (4) 

with U{t) the evolution operator (See Eq. ([3|)). 
The electric-current density of this state is cal- 
culated using the quantum mechanical expression 
j = (Q h/m) 9 [ip*VTp] - ^ 3? bP*M]- After a lengthy 
calculation we arrive at the electric current and separate 
the terms in powers of Aq (hidden in ^). The zero-order 
term does not give new phenomena, and so it is not pre- 
sented. The first-order and second-order terms are 



ji(r,i) 



;f|3{(f)p.o.=--"'/WM5) 



h(r.t) = (^) ^3{/(r)'V/(r)) (6) 

Let us comment on the main characteristics of these two 
currents. ji(r, t) contains the vector momentum ma- 
trix element p„ocO' which reappears in the calculation 
of the current due to the mixing of the two terms of 
i/;(r, t). Thus, as p^oco itself, this component of the 
macroscopic current has a microscopic origin. This cur- 
rent is somewhat analogous to the optical polarization 
in the standard interband optical transitions induced by 
plane waves. Note that this current has no equivalent in 
atomic physics [21| . In Fig [2] we show the current vector 
field for the cases Z = 1 to 4. A clear feature of the cur- 
rent for Z = 1 is that it displays a net circulation around 
the z-axis, which will give rise to a sizable magnetic field. 
For all other values of I the current shows net circula- 
tions around off-centered axes, like in the cases I = 2 to 
4 seen in Fig [21 but no net circulation around the center 
of the beam. Finally, without explicitly showing it here, 
we mention that the current field behaves like a trav- 
elling wave in the z-direction, and therefore, at a given 
time, clockwise and counterclockwise circulations alter- 
nate as one moves along the z-axis. The second-order 
term J2(r,i) comes from just the second term in ?/'(r, i). 
As can be seen in Eq. ([5]), J2(i", stems directly from 
the function /(r), defined after Eq. ([3]), whose spatial 
dependence mimics closely that of laser beam and varies 



significantly only on a macroscopic scale. Thus, we in- 
terpret that this current has a macroscopic origin and is 
analog to the motion induced by twisted light on atomic 
systems. 

For the currents with net circulation around the beam 
axis, e.g. the case / = 1 in Fig. [21 we estimate the mag- 
netic field bz at the center the beam. We adopt a semi- 
classical approach and use Biot-Savart's law with the cur- 
rents from Eq. ([5]) and ([6]) whose azimuthal components 
are 



i20(r,t) 



h 



Ji((7||r||) sin (q^z) t, 
qh 1 



iL^ 



■Mq\\ 



(7) 
(8) 



Then, the total magnetic field is determined by B^ — 
Nex bz- Our semiclassical approach is an admissible pro- 
cedure in perturbation theory since the quantum fluctu- 
ations of the current are of higher order in the vector 
potential and thus smaller than their mean values. 

As an example, we estimate the magnetic field Bz for 
a realistic situation. We choose the following typical ex- 



perimental parameters [IJ, [IJ, |22, |23|. For the laser 
beam: TiiSaphire laser with frequency A ~ 589 nm, rep- 
etition rate F = 100 MHz, power P = 1 /iJ s^^, pulse du- 
ration f = 10 fs, qz ~ 10^ m^^, (j|| = 4 10^ m^^ (spot size 
'^'^^ and 1 = 1. For the semiconductor material 



~ (0.5 ^m)^), and; = 1. 

we choose GaAs parameters: 



10^ 



0.6 nm. 



V = c/3, L = 1.3 //m and \€„ ■ PcOvo\ = 6.4 10^^^ kg m/s. 
Assuming a top-hat laser pulse (which approximates the 
cw-field of Eq. ([21) for long enough pulses) these values 
yield N^x = 2 10^ and (|^| t/h) = 0.95. The resulting z- 
component of the magnetic field in the center of the beam 
is \Biz\ = 1.2 mT and \B2z\ = 0.4 ^tT for the first-order 
and second-order currents, respectively. This estimate 
is reliable since the build-up of the magnetic field takes 
place in a time shorter than the typical relaxation times 
of electrons in the conduction band of semiconductors at 
low excitation density [ISJ]. 

Our theory can be put to test by measuring this 
predicted magnetic field. We suggest the use of time- 
resolved Faraday rotation. A quick estimate of the po- 
larization rotation angle can be obtained with the Vcrdct 
formula 9 = VBzi, where V is the Verdet constant 
and £ is the length traversed by the probe beam. Tak- 
ing V = 10~^rad/(G cm) and i = 10 /xm one obtains 
6 ~ 10~® rad, which is clearly measurable according to 



Ref. [24 1 



In conclusion, in this Letter we have shown for the first 
time the effect that light carrying orbital angular momen- 
tum has on semiconductors. We started by generalizing 
the basic theory of photo-excitation with plane waves to 
the case of Bessel mode beams. We found that the ex- 
citation process generates a superposition of conduction- 
band states with a well-defined angular momentum. This 



superposition state is associated to an electric-current 
density, which acts as the source of a magnetic field. As 
an example, the magnitude of this magnetic field has 
been estimated in a typical experimental situation. The 
important regime of excitonic generation has not been 
addressed here and will be the subject of a future work, 
as well as the inclusion of electron-electron Coulomb in- 
teraction effects. 
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